Scaling properties of granular rheology near the jamming transition 
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Rheological properties of a dense granular material consisting of frictionless spheres are inves- 
tigated. It is found that the shear stress, the pressure, and the kinetic temperature obey critical 
scaling near the jamming transition point, which is considered as a critical point. These scaling laws 
have some peculiar properties in view of conventional critical phenomena because the exponents 
depend on the interparticle force models so that they are not universal. It is also found that these 
scaling laws imply the relation between the exponents that describe the growing correlation length. 



Jamming is a general concept that a disordered sys- 
tem cannot flow under nonzero driving force (e.g., shear 
stress) , which is generally realized at higher densities and 
lower temperatures. Particularly, in a class of athermal 
(i.e., zero temperature) systems such as emulsions, foams, 
and frictionless granular particles, the shear modulus and 
the bulk modulus emerge/ vanish above/below a definite 
density 0, S, i, H . This rigidity transition is now referred 
to as the jamming transition and the transition point is 
referred to as point J. Thanks to massive experiments and 
numerical simulations which indicate the growing corre- 
lation length and characteristic time [I i 0, H H, it 
seems to be confirmed that point J is a critical point, 
although the critical fluctuation is not observed in static 
quantities but found in purely dynamic quantities such 
as the particle displacement in a certain time window. 
Quite interestingly, the critical fluctuation in dynamic 
quantities is also observed in supercooled liquids |ToL [ll| 
so that glass transitions might be related to the criticality 
of point J. Therefore, clarification of the critical nature 
of point J is important in the context of glass physics as 
well as granular mechanics. 

As a result of the criticality of point J, Olsson and 
Teitel find a scaling law for rheology [13], which is of 
the same form as those in conventional critical phenom- 
ena. Such scaling laws provide us a systematic way to 
estimate critical exponents, which are of special impor- 
tance in the context of the universality of critical phe- 
nomena. With regard to the jamming transition, how- 
ever, it is known that some exponents that describe the 
elastic moduli depend on the force model of particles Q 
so that the universality in the sense of conventional crit- 
ical phenomena is questioned. Along the line of thought, 
the model dependence of the critical exponents must be 
further investigated via scaling laws. In this paper, sev- 
eral scaling laws for rheological properties of a granular 
material are investigated by means of molecular dynam- 
ics simulation. It is found that these scaling laws are 
quite different from those for conventional critical phe- 
nomena in the sense that the exponents depend on the 
interparticle force models. We then briefly discuss the 
relation between these exponents and the growing corre- 
lation length. 

We consider a bidispcrsc mixture of frictionless parti- 



cles, the diameters of which are d and 0.7d, respectively. 
The diameter and the position of particle i are denoted 
by Ri and r^, respectively. The force between parti- 
cles i and j is written as Fy = [kF(Sij) + Cn.y • r,,] n. 



where ny = ry 



Tj, and Sij denotes the 



overlap length, (Ri + Rj) — \r t j\. (Note that Sij = if 
(Ri + Rj) < |ry|). Here we test two types of interactions: 



F(z) = z (the linear force model) and F(z) 



,1.5 



(the 



Hertzian force model). Throughout this study, we adopt 
the units in which d = 1, M — 1, and ( = 1. Here M 
denotes the mass of the larger particles, while that of the 
smaller particles is 0.7 3 M. We consider sufficiently hard 
spheres; k — 10 2 for the linear force model and k = 10 4 
for the Hertzian force model. 

In order to ensure uniform shear flow, we adopt 
the SLLOD equations together with the Lees-Edwards 
boundary conditions 14 1. 



q* = — hiQz^y, 

mi 

Pi = ^2 F ii -7Pz,i™y 



(1) 

(2) 



where 7 denotes the shear rate. The system is of con- 
stant volume and consists of 10 3 particles Here we 
investigate several volume densities ranging from 0.5 to 
0.7. We discuss the rheological properties of this system 
with respect to two control parameters: the shear rate 
7 and the volume density <j>. Note that the shear stress 
and the pressure is defined through the virial [3]. 

First we discuss the shear rate dependence of the shear 
stress. As is shown in Fig. [TJ the linear model and the 
Hertzian model show qualitatively the same behaviors, 
albeit quantitatively different. At lower densities, the 
shear stress is proportional to the square of the shear 
rate, which is known as the Bagnold's scaling 15, 16j]. At 
higher density, the shear stress seems to have a nonzero 
value in the zero shear rate limit, which is interpreted as 
the yield stress. From these data, it is found that the 
threshold density <pj (i.e., the jamming density) is some- 
where between 0.645 and 0.648 in the present system. 

As is shown in Fig. [21 these rheological data are col- 
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FIG. 1: The shear rate dependence of the shear stress of 
granular materials subjected to stationary shear flow at sev- 
eral densities, (a) The model in which the particles interact 
via linear repulsive force. (The linear force model), (b) The 
model in which the particles interact via Hertzian force. (The 
Hertzian force model). 



FIG. 2: Collapse of the rheological data show in Fig. [T] using 
the scaling law, Eq. (|3]). The solid lines are proportional to 
j 2 , which represent the Bagnold's scaling, (a) The scaling 
law Eq. (|3} for the linear force model, (b) The scaling law 
Eq. ((3]) for the Hertzian force model. 



TABLE I: The exponents estimated from scaling laws, Eqs. 
©, ©, and @. 



Observable 


Shear stress Pressure Temperature 


Exponent 


ys> yt/vi v'^lv'-, ^* x#/a; 7 


Linear force 
Hertzian force 


1.2(1) 1.9(1) 1.2(1) 2.1(2) 2.5(2) 1.9(2) 
1.8(1) 2.4(1) 1.8(1) 2.4 (1) 3.2(2) 2.4(2) 



lapsed using the following relation. 



s= m y *s± 



7 



|$|b*/»7 



(3) 



where S± (■) is the scaling function and $ denotes cpj — <fi. 
The jamming density <pj is set to be 0.646, which we use 



throughout the present numerical analysis [17[ • Note that 
this scaling law, Eq. ([3]), is of the same form as those in 
conventional critical phenomena. One may thus expect 
that the jamming transition is one of the conventional 
critical phenomena and that the exponents y$ and y 7 , 
together with the scaling dimension of the shear stress, 
classify the jamming transition into a universality class. 
However, this idea is ruled out because these exponents 
depend on the force model. We obtain j/$ = 1.2 ± 0.1 
and y<s>/y 7 = 1.9 ± 0.1 for the linear force model, while 
yg, = 1.8 ± 0.1 and y<s>/y 7 — 2.4 ± 0.1 for the Hertzian 
force model. The exponent y 1 is thus roughly estimated 
as 0.63 and 0.75 for the linear model and for the Hertzian 
model, respectively. Recalling that the elastic moduli of 
granular media depends on the interaction force [3j, the 
model dependence of ?/$ is rather obvious in view of Eq. 
(j4|). Furthermore, note that the other exponent y 1 also 
depends on the force model. These exponents are listed 
in Table. |U 



The scaling law Eq. (J3]) implies two nontrivial rheo- 
logical properties. First, it describes power-law rheology 
at $ = 0, where the right hand side of Eq. ([3]) does not 
depend on <£>. This leads to S oc 7 y T. Note that the 
exponent y 7 ~ 0.63 for the linear force model is directly 
observed in the previous studies [3, Qji]. Second, Eq. 
([3]) also describes the density dependence of the dynamic 
yield stress. Inserting 7 = into Eq. ([3]), we get 



S($,0) = S_(0)|$| 



(4) 



The dynamic yield stress is thus proportional to |<I>| ^/ ' , ', 
provided that <S_(0) ^ 0. 

Then we discuss the behavior of the pressure (i.e., the 
normal stress). One of the essential features of gran- 
ular rheology is that the pressure also depends on the 
shear rate [18| , because granular systems have no intrin- 
sic temperature and the extent of the random motion of 
particles increases with the shear rate. Remarkably, they 
also obey a scaling law; 



7 



(•5) 



where V±(-) is the scaling function. The exponents are 
estimated as = 1.2 ± 0.1 and y$/y 7 = 2.1 ± 0.2 for 
the linear force model, and y'^ = 1.8 ± 0.1 and y'^/y'^ = 
2.4 ± 0.1 for the Hertzian force model. The exponent 
y' 7 for each model is thus roughly estimated as 0.57 and 
0.72, respectively. Note that the exponents again depend 
on the force model. 

This scaling law, Eq. ([5]) , again implies two important 
rheological behaviors. First, inserting 7 = into Eq. 
([5|) leads to the density dependence of the pressure in 
the 7 — > limit; P(<&,0) oc l^*. We remark that this 
pressure is the dynamic yield pressure and thus should 
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FIG. 3: The shear rate dependence of the pressure (i.e., the 
normal stress P zz ) of granular materials subjected to station- 
ary shear flow at several densities, (a) The linear force model, 
(b) The Hertzian force model 



10" 



± 10 
ft 

10" 
10" 
10"' 

io" : 




FIG. 5: Collapse of the shear rate dependence of the kinetic 
temperature using the scaling law, Eq. (B). The kinetic tem- 
perature is defined as T — (iriiSvf) /3. In the both panels, 
the solid lines are proportional to j 2 , which represent the 
Bagnold's scaling. The dotted lines are proportional to 7, 
which indicate that the kinetic temperature is proportional 
to 7 above the jamming density, (a) The linear force model, 
(b) The Hertzian force model. 



scaling properties of the kinetic temperature defined by 
T = (mjtfvf )/3. Note that this kinetic temperature is dif- 
ferent from several effective temperatures that have been 
tested in sheared athermal systems, which take a nonzero 
value in the static limit [20| . These effective temperatures 
are dominated by the static limit value and hence insen- 
sitive to the shear rate, while the kinetic temperature is 
much more sensitive to the dynamics. The scaling law 
for the kinetic temperature reads 



FIG. 4: Collapse of the rheological data shown in Fig. |3]using 
the scaling law, Eq. ((5j). The dashed lines are proportional 
to j 2 , which represent the Bagnold's scaling, (a) The scaling 
law for the linear force model, (b) The scaling law for the 
Hertzian force model. 



not be identified with the pressure at zero shear strain, 
which is proportional to |<1>| 10 for the linear force model 
and I $ 1 1,5 for the Hertzian force model 0, Sj- Second, 
inserting $ = into Eq. §5§ leads to the power-law 
behavior at the jamming density; P(0, 7) oc y y i. Note 
that y' ~ 0.57 for the linear force model is consistent 
with the previous study [H, [l9| . 

As is discussed above, the exponents in the scaling 
laws depend on a force model. This is partially under- 
stood if one recalls that y$ and y'^ describe the dynamic 
yield stress and the dynamic yield pressure, respectively. 
These quantities are essentially dominated by the elastic 
moduli of the system, the exponents for which do de- 
pend on the force model 0. Along the line of thought, 
dynamical quantities that vanish in the static limit may 
have a universal scaling behavior which does not depend 
on the force model. To test this idea, we discuss the 



T($, 7 ) = |$|**T ± 



7 



(6) 



where 7±(-) are the scaling functions. We stress that 
the exponents again depend on the force model: x$ = 
2.5 ±0.2 and x^/xj = 1.9 ±0.2 for the linear force model 
and x$ = 3.2±0.2 and x<& /x 1 = 2.4±0.2 for the Hertzian 
force model. Therefore, the exponents for the kinetic 
temperature, which is a purely dynamic quantity, do de- 
pend on the force model. Therefore we can conclude that 
the exponents in the rheological scaling laws are not uni- 
versal in contrast to conventional critical phenomena. 

The scaling laws that are obtained in the present sim- 
ulation, Eqs. ([3]), ([5]), and ©, involve several exponents, 
the physical meaning of which is still not clear. In the rest 
of this paper, we discuss two important points regarding 
their physical meaning and their universality focusing the 
criticality of point J. First, we discuss the model depen- 
dence of the exponents in the scaling law for rheology, 
Eq. ([3]), by comparing it with the scaling law obtained 
by Olsson and Teitel (OT) \v^. 



S 

W 



(7) 



4 



where OT obtain A = 1.2 ± 0.1 and (3 = 1.7 ± 0.2 for a 
two dimensional system with the linear force model. In 
order to see the correspondence between Eqs. fl3J) and 
(|7|), it is more convenient to rewrite Eq. as 



1*1 



9± 



7 



|$|A+/3 



(8) 



Now it is apparent that A is identical to j/$ in Eq. |3|), 
and actually j/$ is 1.2 ± 0.1 for the linear force model 
in the present three dimensional system. Therefore, the 
exponent A (or ?/$) is independent of dimensionality, but 
depends on the force model. This is plausible when we 
recall the behavior of the shear modulus with respect to 
the density is independent of dimensionality but depends 
on the force model We also remark that the other 
exponent (3 is not the same for the two cases: j3 ~ 1.7 
[l2| and ~ 0.7 here. This discrepancy results from par- 
ticle dynamics; OT adopt overdamp (massless) dynamics 
whereas the inertia comes into play in the present sim- 
ulation. We indeed obtain j3 ~ 1.7 for a three dimen- 
sional system with overdamp dynamics, which is to be 
presented elsewhere. 

The second point we wish to discuss is the relation 
between the scaling laws and the exponents that describe 
the growing correlation length. To see this, we rewrite 
Eq. ([3|) in the following form. 



5(*,7) = A- 1 5 , (A 1 / 2 "-$,A 1/ ^ 7 ), 



(9) 



where A is an arbitrary nonnegative number. Here we as- 
sume that the correlation length £ diverges as the system 
approaches point J; i.e., $ — > and 7^0 (or T — > 0). 
Then, following the idea of conventional critical phenom- 
ena, Eq. © implies that £ - S- ys , £ ~ |$|~^ s , and 



ViVS 



, where 1 jys is the scaling dimension of the 



shear stress. These relations indicate that the exponents 
that involve the algebraic divergence of the correlation 
length are not independent. Here we only discuss the 
exponents for the linear force model, because few results 
are known as to the growing correlation length in the 
Hertzian force model. Using £ cx | — ^ | 0.25 ^at been 
confirmed in a sheared granular material (l9j . we obtain 
ys ~ 0.4 so that £ ~ |<I>|~ - 5 . Note that this exponent 0.5 
agrees with the result of Wyart et al. 2l| , which involves 
the mean free path of force chains. We remark that the 
same result can be obtained using Eqs. ^ and ^ in- 
stead of Eq. ([3]), in which case the scaling dimensions 
of the pressure and the kinetic temperature are approxi- 
mately 1/0.4 and 1/0.2. respectively. 

To summarize, we find the scaling laws that describe 
the shear stress, the pressure, and the kinetic tempera- 
ture near the jamming transition point. These scaling 
laws are related to the exponents that describe the grow- 
ing correlation length, and hence advocate the critical- 



ity of point J. However, the 'critical' exponents is found 
to depend on the force model, which indicates the pe- 
culiarity of point J in contrast to conventional critical 
phenomena. 

The author gratefully acknowledges helpful discussions 
with Michio Otsuki, Hisao Hayakawa, and Namiko Mi- 
tarai. 
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